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ABSTRACT 

We present the first study of the dynamical evolution of a star cluster that combines a significant 
population of primordial binaries with the presence of a central black hole. We use direct N-body 
simulations, with a black hole mass of about one percent of the total mass of the cluster. 

The evolution of the binary population is strongly influenced by the presence of the black hole, 
which gives the cluster a large core with a central density cusp. Starting from a Plummer profile, 
we first encounter a phase, that last approximately 10 half-mass relaxation times, in which binaries 
arc disrupted faster compared to analogous simulations without a black hole. Subsequently, however, 
binary disruption slows down significantly, due to the large core size. 

The dynamical interplay between the primordial binaries and the black hole thus introduces new 
features with respect to the scenarios investigated so far, where the influence of the black hole and 
of the binaries have been considered separately. Specifically, the pattern of binary destruction by an 
intermediate-mass black hole may leave a fingerprint that could be detected observationally. 
Subject headings: stellar dynamics — globular clusters: general — methods: n-body simulations — 
binaries: general 
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1. INTRODUCTION 

Over the last few years some tantalizing, but yet far 
from conclusive, evidence has been accumulating in sup- 
port of the idea that some star clusters could harbor 
a central black hole (BH) with a mass of the order 
of 10 3 M Q or more. Detection of such an intermediate 
mas s black hole (IMBH ) has been claimed for M15 and 
Gl llCerssen et aljl200l ICehWdt et al .ll2fM Elf!! . al- 
though accurate dynamical models of Ml 5 and Gl can 
be obta ined without a central BH <)Baumgardt et af] 
l2003albi ). Interestingly, the visual appearance of glob- 
ulars containing an IMBH is not that of a so-called core- 
collapsed clu ster, but rather that of a cluster with a still 
sizable core l|Baumgardt et aLl |2005). 

IMBHs present a high theoretical and observational in- 
terest as these could be potential ultra-luminous X-ray 
sources and even emit gravitational waves, detectable by 
the next generation of gravitational wave detectors, as a 
result of close interactions with stars. However, despite 
this interest and the fact that theoretical studies of BHs 
in stellar systems started more than 30 years ago (e.g., 
see lPeeblesllT97l iBahcall fc WolfllT97^ . detailed direct 
N-body simulations to study the dynamics of an ideal- 
ized model with single st ars and a central BH have b een 
performed only recently ijBaumgardt et aLll2004albt see 
also the studies on the for mation of IMBHs by runawa y 
mergers of massive stars bv lPortegies Zwart et al.l [2004'i. 

One ingredient that complicates N-body simulations 
of globular clusters is the presence of primordial bi- 
naries. Often, these are neglected in large simula- 
tions despite the increasing observational evidence that 
many star s in a globular clus t er ha v e a companion 
ftTiit et, alJ ITflfll lAlhrow et a.U HQ IBellazzini et a,D 
12002 iZhao fc Bailvnl |2005|) . This frequent neglect 
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is due to the dramatic increase in computational re- 
sources required in a simulation where the local dy- 
namical timescale may be many orders of magnitude 
smaller than the global relaxation timescale (hard bina- 
ries have an orbital period of a few hours, while the half- 
mass relaxation time can be up to a few billion years). 
The study of the dynamics in the presence of primor- 
dial binaries has been mainly limited to Fokker-Planck 
or Monte Carlo a pproaches l|Giersz fc Spurzeml 1200(1 
iFregeau et al.l2003T) and to direct simula tions with rather 
modes t particle number s , from N sa 10 3 ([M cMillan et^dj 
Il990t IMcMillan fc Hutl fl99l IHeggie fc Aarsethl 119921) 
to recent higher resolution s imulations, with N up 
to 16384 ijHeggie et alJ 12005(1 . Some realistic sim- 
ulations, including primordial bi naries are available 
(jPortegies Zwart fc McMillanl 12004]) . but these are lim- 
ited only to the first stage of the life of young dense 
clusters. In the case of open clusters, M67 has been 
modeled in a 36,000-body sim ulation running for several 
Gigayears HHurlev et al.ll2005h . 

The presence of either an IMBH or a significant popu- 
lation of primordial binaries leads to an early release of 
abundant energy, inhibiting the development of a deep 
core collapse and hence of the onset of gravothermal os- 
cillations. An IMBH can generate energy by swallowing 
or tightly binding stars deep in its potential well [note 
also that energy can be generated through encounters of 
stars in the density cusp that is formed around the BH] , 
while primordial binaries can generate energy by rapidly 
increasing binary binding energy through three and four 
body encounters. Energy thus generated in the core of 
the system fuels the expansion of the half-mass radius, 
leading to a self-similar expansion of the entire system 
<|Henonll19fl5D . 

In this work we present the first direct simulations of 
the evolution of globular clusters with both a significant 
fraction of primordial binaries (10 %) as well as an IMBH 
with a mass of 1.4 — 2.5 % of the total mass of the sys- 
tem. We address the following questions. Under the 
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combined effect of the BH and of the binaries, what is 
the equilibrium size of the cluster core? How is the bi- 
nary population affected by the presence of the central 
IMBH? Which physical processes dominate in the core? 
In principle two competing effects are possible: the BH 
may either enhance the disruption rate of binaries both 
indirectly, due to the cre ation of a de nsity cusp, and di- 
rectly, by tidal stripping (Pfahl 2005), or it may reduce 
the probability of interactions between binaries and sin- 
gles, by pr oducing a low stellar de nsity in a relatively 
large core l|Baumgardt et al1l2004ajL Which process is 
dominant is of fundamental importance. These questions 
are addressed in the next sections. 

2. NUMERICAL SIMULATIONS: SETUP 

The simulations presented in this paper have bee n per- 
formed using the NBODY-6 code ijAarsethl 12003). that 
has been modified with the kind help of Dr. Aarseth 
to ensure a more efficient and accurate treatment of the 
dynamics around the BH (the relative energy error at 
the end of our simulations is below 0.5 %). We used 
a total of N to t — 9011 equal mass stars (each of mass 
m), arranged in 819 binaries and 7373 singles, plus a 
BH, introduced as a massive star, with mass (tush) in 
the range 1.4 — 2.5 % of the total mass of the system. 
The system is isolated and the initial binding energy dis- 
tribution for the binaries is in the range 3 to 400 kT 
(where 3/2 kT is the average kinetic energy of a star 
of the system), f lat in logarithmic scale as in the runs 
bv iHeggie et alJ l|2005|) . that we use for comparison, as 
here we employ the same initial conditions except for 
the presence of the BH. The initial mass distribution is 
that of a Plummer model. To initialize the simulation 
in a situation of approximate dynamical equilibrium in 
the presence of the BH we generate a Plummer model 
made of single stars only, we then scale the velocities of 
the particles to reach virial equilibrium and let the sys- 
tem evolve for 5 half-mass crossing times. Some stars are 
then selected at random and binaries are added to the 
simulations, with the standard initialization provided by 
NBODY-6. The evolution of the system is then followed 
up to t ps 25 t r h{0)- For t r h(0) we mean t r h com puted at 
t = 0. Here t r h is the half-mass relaxation time l)Spitzerl 

IT987T) . defined as t rh = 0. 138 Nr^ 2 /(VMG In (0.11AT)), 
with Th being the half-mass radius, M the total mass of 
the system and N the number of centers of mass (i.e. 
N = 8192 in our case); in NBODY units t rh (0) w 112 
for our initial conditions. 

3. GLOBAL EVOLUTION 

The large scale structure evolution of the cluster is 
dominated by the heating related to the presence of the 
BH. In fact only the inner regions of the system expe- 
rience a mild collapse on a timescale that, depending 
on the mass of the BH, is of the order of a few t r h- 
Inside the c ore radius r ri defined a s the density aver- 
aged radius l)Casertano fc Hutlll985fL a cusp in the den- 
sity profile is formed within t he sphere of influence ry 
(with Ti w 15 r c rriBH/M ', see iBaumgardt et al.l l2004a) 
of the BH, with a profile proportiona l to ~ 1/r 1 - 7 and 
thus sim ilar to the l/r 1 - 75 measured bv lBaumgardt et alJ 
(2004a). For tush = 0.014M the influence radius is ap- 
proximately 0.2r c . By definition, the stellar mass within 



TABLE 1 
Summary of runs with N = 8192 



m bh /M f (r c )cc/r c0 r c /r h r hf /r h0 

0.014 0.1 0.75 0.29 2.46 

0.025 0.1 0.80 0.31 2.83 

0.014 0.70 0.27 1.92 

0.1 0.25 0.09 2.25 

Note. — In the first column we report the BH mass, in the 
second the fraction / of primordial binaries, in the third the core 
radius at the end of the initial core contraction phase (r c ) cc in units 
of the initial core radius r c o, the fourth entry is the core to half 
mass radius ratio during the self-similar expansion of the system, 
while the last entry {r^f) is the value of the half mass radius at 
t = 24 t r h (0) in units of the initial half mass radius. 



this radius is comparable to that of the BH, and thus 
around one percent of the total mass of the cluster. 

In case of ttibh = 0.014M, the core radius r c is re- 
duced in ps 4 t r ho from the initial value of 0.4 to 0.3 (in 
NBODY units). This is to be compared with a value 
of 0.1 reached without the BH; when a BH (with the 
same mass) but no primordial binaries are present r c goes 
down to ps 0.28. 

After the first mild contraction all Lagrangian radii 
start to expand steadily and a self-similar regime sets 
in, with the half-mass radius growing in proportion to ps 
t 2 / 3 , in agreem ent with the theoretical argument given by 
iHenonl l)1965fl . The half-mass radius is thus marginally 
bigger (by « 10 % at 24 t r ho) due to the presence of 
the BH; without a BH the half-mass expansion starts 
only after core collapse, which takes ps 10 t r ho in that 
case. Conversely, if a BH but no binaries are present, the 
expansion rate of is reduced by ps 20 %. A summary of 
the properties of our simulations is reported in Tabled 

4. PROPERTIES OF THE BINARY POPULATION 

As the presence of the BH dominates the global dynam- 
ics of the cluster, the evolution of the binaries presents 
some remarkable differences from the scenario where no 
central BH is present. We can distinguish two phases. 

The presence of a BH significantly accelerates the rate 
at which binaries are disrupted in the first few half-mass 
relaxation times (see Fig.^). The reason is that binaries 
that pass near the BH can be quickly destroyed. The 
disruption rate is approximately constant during the first 
5 t r ho, while without a BH disruption takes a while to 
get underway. The initial disruption rate depends on 
the mass of the BH: a more massive BH starts burning 
binaries at a higher rate, as can be seen from Fig. ^ 

For our simulation with msn = 0.014M we observe a 
binary depletion rate of ps 0.15 binaries per initial half- 
mass crossing time, during the first 5 t r ho (80 binaries in 
total) . Around half of the binaries are disrupted within 
n . Interestingly only two of them happen to be so close 
to the BH tidal radius r± that the disruption may be con- 
sidered a direct consequence of the tidal stripping force 
exhorted by the BH. Here r t — (tyibh I Tfi)^'^a ps 5a, 
where a is the semi-major axis of the binary. Wit h this 
respec t, the theoretical model recently proposed bv lPfahj 
(2005) applied to our simulation would give (from his 
Eqs. 11-12) a tidal disruption rate of ps 2 • 10~ 2 binaries 
per initial half-mass crossing time, so that we observe 
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Fig. 1. — Time-dependence of the number of single (upper set of 
lines) and binary stars (lower set of lines) expressed as a fraction of 
the initial values. The data refer to simulations with N = 8192 and 
10% primordial binaries. There are three choices for BH mass: no 
BH (dotted line), a central BH with msH = 0.014M (solid line), 
and ragu = 0.025A/ (dashed line). The unit of time is the initial 
half-mass relaxation time. 



less tidally stripped binaries than expected (« 10). The 
disruption of binaries by the BH is thus mainly due to 
an indirect effect: binaries venturing close to the BH, 
where the density is higher, are more likely to interact 
with a single star or with another binary. In fact in our 
simulations we often observe the presence of hierarchi- 
cal systems within , with the BH playing the role of a 
perturber. 

After the initial transient phase a self-similar expan- 
sion sets in, where the average core density is much less 
(approximately by a factor 10) than it would be without 
the presence of a central compact object. In this second 
phase we observe (see Fig. ^) a reduced rate of binary 
disruption (as this rate is proportional to the s quare of 
the density, e.g. see iVesperini fc Chernofllll994j) . so that 
the difference between a simulation with and without a 
BH becomes remarkable. The turning point is around 
10 t r ho, when the number of surviving binaries for a sim- 
ulation with a central object becomes greater than in 
the absence of a BH. Our Fig. ^ has been given in units 
of the initial half-mass relaxation time, but the picture 
remains qualitatively the same even if we consider a co- 
moving time coordinate to take into account the differ- 
ences between the half-mass radii of a simulation with 
and without a BH. 

Interestingly the spatial distribution of binaries is also 
different from the case without a BH: as can be seen from 
FigEl the number of binaries within the 0.05 Lagrangian 
radius is much less for the simulation with a BH. This is 
probably due to the disruption of binaries that approach 
close to the BH. In the presence of a BH, binaries tend 
to be more concentrated between the 0.05 Lagrangian 
radius and the half-mass radius, while in the absence of 
a BH, binaries can sink deeper into the central region of 
the cluster. As a result, the ratio of binaries to singles 
in the core of the system is between two and three times 
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Fig. 2. — Distribution of binding energies for binaries (upper 
panel) at time t = 24 t r h(0) for a simulation with N = 8192, 
10 % primordial binaries and mgu = 0.014M (solid line), com- 
pared to a similar simulation without the BH (dotted line). In the 
lower two panels we report the energy-radius distribution of bina- 
ries for the same simulations. The radii have been normalized to 
the Lagrangian radius containing 5 % of the mass in the simula- 
tion. The dashed line gives the position of the half-mass radius, 
at t = 24 t r i l ( l . The relatively large survival at t = 24 t r ^o of bi- 
naries within the half-mass radius with binding energy Ei, < 20kT 
is clearly visible in the bottom panel, and is also reflected in the 
higher values for the solid line at the left-hand side of the upper 
panel. 



less than in absence of the BH. 

If we compare the evolution of the binding energy dis- 
tribution of the binaries in a run with and without a 
central BH, we can see that the net effect of the BH at 
later times is to somewhat enhance the survival prob- 
ability of binaries with binding energies of a few kTs 
(Ef, < 16 kT), especially for binaries within the half-mass 
radius (see Fig. [2J . This feature is especially apparent in 
a two-dimensional plot of the distribution of binaries as 
a function of binding energy and radial distance from the 
cluster center, depicted in the bottom panels of Fig. [21 
The rather strong correlation betwee n radius and bind- 
ing energy that is observ ed (e.g., see iGiersz fc Spurzeml 
2000; Hcg gie et alJ 120051) in systems without a central 
compact object is thus almost absent. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have presented for the first time the 
results of direct N-body simulations of the evolution of a 
globular cluster with a population of primordial binaries 
that harbors a central IMBH, with a mass of the order 
of 1 % of the total mass of the system. In order to be- 
gin to analyze the basic evolutionary processes we have 
restricted our simulations to isolated systems of stars of 
equal mass, while for the time being neglecting stellar 
evolution and physical collisions. 

The environment around the BH turns out to be of 
great interest from a dynamical point of view. It repre- 
sents a laboratory where complex interactions between 
hierarchical systems take place. Around the BH we ob- 
serve usually one or more stars tightly bound to it. With 
a significant population of primordial binaries, it fre- 
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quently happens that a binary approaches the center of 
the system. Interactions between these subsystems are 
important not only from a dynamical point of view, but 
also because they form a factory to produce exotic ob- 
jects, such as tight (X-ray emitting) binaries and high 
velocities escapers, which we observe in high numbers in 
our simulations. In our run with tubh = 0.014M we find 
that a fraction of 3 • 1CP 3 of the stars of the cluster leave 
the system in a relaxation time and that « 1/3 of them 
have escape velocities that exceed by at least five times 
the core velocity dispersion. Of these 54 were binaries; 
in addition one triple was ejected from the system. 

From a computational point of view the complex in- 
teractions between binaries and the BH have proved to 
be very challenging for current state of the art N-body 
algorithms. Clearly, new methods have to be developed 
to provide both more accuracy and higher efficiency. 

We have shown that a BH with a mass of around 1 % of 
the cluster mass can produce most of the energy required 
to fuel the expansion of the globular cluster, so that the 
evolution of the large scale structure of a globular cluster 
with primordial binaries and a central BH is remarkably 
similar to that attained in the absence of binaries (see 
iBaumgardt et aLl l2004a'). Even so, the evolution of the 
population of primordial binaries is strongly influenced 
by the BH, so that in principle indirect evidence for the 
presence of the BH can be obtained by observing the 
distribution of binaries. 

Observational detection of an unexpectedly large num- 
ber of only moderately hard binaries in the core could be 
used as arguments for the presence of a IMBH in old star 
clusters, i.e. in those systems whose age would imply a 
depletion in the central region of binaries with binding 
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energy below w 15 kT. These binaries are able to survive 
in significant numbers when an IMBH is present. 

Clearly, a quantitative comparison with observations 
would require us to proceed beyond the simple models 
presented here. To start with, we would need to include 
a realistic number of stars to avoid possible biases in- 
troduced by the scaling of the ratio of the BH mass to 
that of single stars. In addition, the introduction of a 
realistic initial mass function is likely to modify the con- 
centration of binaries in the center of the system due 
to mass segregation, and stellar evolution will also influ- 
ence the distribution of binary binding energy. However, 
the main effects presented in this letter are likely to be 
present, at least in qualitative ways, in more detailed 
realistic simulations. 
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